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X-Ray Photon Correlation Spectroscopy Reveals Intermittent Aging Dynamics in a 

Metallic Glass 
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We use coherent X-rays to probe the aging dynamics of a metallic glass directly on the atomic level. 

Contrary to the common assumption of a steady slowing down of the dynamics usually observed 
in macroscopic studies, we show that the structural relaxation processes underlying aging in this 
metallic glass are intermittent and highly heterogeneous at the atomic scale. Moreover, physical 
aging is triggered by cooperative atomic rearrangements, driven by the relaxation of internal stresses. 

The rich diversity of this behavior reflects a complex energy landscape, giving rise to a unique type 
of glassy-state dynamics. 


Physical aging is not only of scientific interest mia 
but also of great practical importance, as the perfor¬ 
mance stability of many technologically relevant amor¬ 
phous materials such as oxide glasses, polymers and 
metallic glasses (MGs) depends on how aging affects their 
properties during the expected service life. Eor MGs in 
particular, aging can result in severe embrittlement [ 3 ] 
and bring about profound changes to fracture and fa¬ 
tigue properties [ 4 ]. Therefore, the ability to understand 
and characterize the micromechanisms of aging in MGs 
will play a central role in ensuring their success in a broad 
range of industrial and commercial applications. 

A key to understanding aging on the microscopic level 
is measuring the time scale on which the system rear¬ 
ranges its internal structure. These structural rearrange¬ 
ments are directly related to how fast microscopic den¬ 
sity fluctuations in the system decay [ 5 ]. The physical 
quantity describing this is the density correlation func¬ 
tion, which can be measured over a wide range of time 
scales using a variety of complimentary experimental 
techniques [ 6 ]. Additionally, it can be formally defined 
within the framework of statistical mechanics, allowing 
for detailed comparisons between experiment, theory [ 7 ] 
and simulation [ 8 ] . Access to the density correlation func¬ 
tion is therefore fundamental for describing the dynami¬ 
cal behavior of glasses and achieving a consistent picture 
of structural relaxation on microscopic length scales. 

X-ray photon correlation spectroscopy (XPCS) has re¬ 
cently emerged as a novel technique for studying the mi¬ 
croscopic dynamics of condensed matter [H IQUTT]. In 
XPGS, coherent beams of X-rays are scattered from a 
sample and the resulting fluctuations in intensity corre¬ 
lated over time. This enables measurement of the den¬ 
sity correlation function on the atomic scale. Previous 
XPGS investigations of hyper-quenched Mg65Gu25Yio 
and ZreeNiss MGs m m have revealed the presence 
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FIG. I. Structural relaxation times measured with DMA 
(filled circles) and XPCS during slow cooling from the super¬ 
cooled liquid (filled squares) and subsequent re-heating from 
the partially annealed glass (open squares). The solid line is a 
fit of the equilibrium data to the VFT equation. The dashed 
line is an Arrhenius fit to the XPCS cooling data in the glass. 

of a unique stress-dominated dynamics characterized by 
subsequent aging regimes: a fast exponential growth for 
short waiting times, followed by a long, almost station¬ 
ary regime. The latter has been observed also in network 
glasses m and is contrary to the steady slowing down 
of the dynamics during aging observed in macroscopic 
studies. 

Here, we investigate the atomic motion in a 
Pd43Cu27NiioP2o MG ribbon using XPGS and show that 
this system ages in a highly heterogeneous manner, con¬ 
sisting of quiescent periods of stationary dynamics in¬ 
terspersed with cooperative, avalanche-like structural re¬ 
arrangements. Our findings suggest the existence of 
a complex dynamical scenario occurring at the atomic 
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level and are consistent with a stress-dominated dynam¬ 
ics that controls aging. The XPCS measurements were 
performed at the beamline IDIO at ESRF, France. A 
detailed overview of the experimental set-up and data 
treatment procedure can be found in the Supplementary 
Material m- 

The structural relaxation times r of Pd 43 Cu 27 NiioP 2 o 
in the ultra-viscous supercooled liquid above Tg were 
determined from dynamic mechanical analysis (DMA) 
measurements [14] employed at different test frequen¬ 
cies and are given by the filled circles in Fig. We 
model the temperature dependence of r by fitting the em¬ 
pirical Vogel-Fulcher-Tammann (VFT) equation r(T) = 
To exp[D*To/(T — To)] to the equilibrium data. Here D* 
is a measure of an inherent liquid-state property, the so- 
called fragility [15]. This reflects the degree of depar¬ 
ture of the liquid dynamics from the Arrhenius equation, 
with smaller values of T* corresponding to more fragile 
liquids. We find T* = 10.3 ± 0.4, in agreement with 
previously reported rheological data [T 6 ] . 

The square symbols in Fig. [^are the values of r mea¬ 
sured using XPCS. Below Tg = 563 K, the transition into 
the out-of-equilibrium glass is apparent. In the glass r 
displays a much weaker temperature dependence than in 
the supercooled liquid which is well described by an Ar¬ 
rhenius equation with an activation energy of 178.9 ± 0.4 
kJ/g-atom. The open squares represent r measured dur¬ 
ing subsequent reheating of the partially aged glass from 
533 K. Relaxation in the glass depends on the previous 
thermal history, reflected here by the higher values of r 
as a result of the sample aging. 

The time evolution of the microscopic dynamics dur¬ 
ing aging is directly captured in XPCS by the intensity 
two-time correlation function G(g,^ 1 ,^ 2 ) (TTCF). This 
quantity reflects the statistical similarity between aver¬ 
age configurations measured at times ti and ^2 and is 
thus the time-resolved version of the standard intensity 
correlation function g 2 {Q^t) = (G) [ 6 ]. The TTCF mea¬ 
sured upon cooling from the supercooled liquid during an 
isothermal hold of 2.2 hr at 548 K {Tg —15 K) is shown in 
Fig. Ha). Along the center diagonal, which corresponds 
to the elapsed time of the measurement, the values of 
the TTCF are largest and decrease as the difference be¬ 
tween ti and t 2 grows. The width of the diagonal contour 
is thus directly proportional to r, i.e. the characteristic 
time scale on which any given atomic configuration no 
longer corresponds to that measured at a later time. 

The normalized ^ 2(^0 5 ^) function, extracted from the 
TTCF after progressively longer waiting times is 
shown in Fig. Hb). Fach correlation function was fit¬ 
ted using a Kohlrausch-Williams-Watts (KWW) expres¬ 
sion g2{Qo^t) = 1 + cexp[— 2 (t/r)^], where r is the relax¬ 
ation time, (3 is the shape parameter m and c = 7/2 
is the product between the experimental contrast 7 and 
the square of the Debye-Waller factor fg. The temporal 
dependence of r and [3 is displayed in the inset. 


(a) 

1 1.014 
1.012 
1.010 
1.008 
1.006 
J 1.004 

1 1.002 
1.000 



(b) 



FIG. 2. (a) Two-time correlation function (TTCF) at T = 
548 K, measured in the glass upon step-cooling from the su¬ 
percooled liquid. The width of the diagonal intensity profile 
is proportional to r. (b) Normalized ^^ 2(^0 = 2.8 A“^,t) ex¬ 
tracted from the TTCF after different waiting times Solid 
lines through the data are best-fits to the KWW equation. 
The values of r and /3 are given in the inset as a function of 
tw A stationary dynamics persists up to 4.6 x 10^ s, after 
which the dynamics slows down as sample ages. 


Generally, aging is expected to begin as soon as the 
sample falls out of equilibrium and r becomes compara¬ 
ble to the experimental observation time [18]. Surpris¬ 
ingly, we observe the onset of aging here only after a long 
stationary regime of 4.6 x 10^ s [Fig. I^b)] . A similar sta¬ 
tionary dynamics has been reported in previous XPCS 
studies of as-quenched MGs laill]. In those works, how¬ 
ever, stationary dynamics was only observed following an 
initial, fast aging regime, while the measurements pre¬ 
sented in Fig. I^show the opposite trend. Nonetheless, 
in the potential energy landscape (PEL) approach [19], 
both stationary dynamics suggest the presence of deep 
energy basins, in which the system becomes trapped and 
explores local minima separated by similar energy bar¬ 
riers. During aging, the system eventually escapes the 
basin and begins to explore deeper energy minima. The 
discrepancy in the occurrence of our stationary dynamics 
and that reported in Refs. mm is likely the simple con¬ 
sequence that our glass is studied upon slow cooling from 
the supercooled liquid, whereas those studies had mea- 













3 



FIG. 3. (a) Measured TTCF at T = 538 K after step-cooling 
from 548 K. The overall broadening of the intensity profile re¬ 
flects a slowing down of the microscopic dynamics correspond¬ 
ing to physical aging, (b) Distribution of r during stationary 
dynamics (T = 548 K) and heterogeneous aging dynamics (T 
= 538 K). 

sured hyper-quenched ribbons trapped already at much 
higher levels on the PEL. 

These observations are consistent with recent molec¬ 
ular dynamics (MD) simulations of a Cu 56 Zr 44 MG by 
Fan et a/., which show that the microscopic relaxation 
dynamics is significantly more spatially localized in a re¬ 
laxed or slowly cooled glass, than in a fast quenched ma¬ 
terial [20] . 

The peculiar dynamical scenario observed at 548 K be¬ 
comes even more dramatic at lower temperatures, where 
a highly heterogeneous relaxation dominates the ob¬ 
served aging behavior. Figurej^a) shows the TTCF mea¬ 
sured after step-cooling from 548 to 538 K and holding 
isothermally for 3.9 hr. The broadening trend of the in¬ 
tensity profile with increasing clearly shows that the 
dynamics is slowing down significantly during annealing. 
A remarkable feature, however, is that the aging pro¬ 
cess cannot be solely described by a more or less steady 
slowing down of the dynamics, as observed in previous 
XPCS investigations [3 US- Instead, rapid decorrela¬ 
tion events - characterized by a sudden narrowing of the 
TTCF intensity profile - appear intermittently as the 
measurement progresses. A closer inspection of ^ 2(^0 5 ^) 
as a function of also reveals the presence of multi¬ 


ple stationary regimes m- The histograms in Fig. |^b) 
show the distribution of r at both 548 and 538 K, ob¬ 
tained from fitting the KWW-equation to the extracted 
92 {Qo^t) [14]. The stationary dynamics at 548 K (Fig. 
is shown by the unimodal behavior of r around 1.5 x 1^ 
s, while the heterogeneous aging at 538 K exhibits a bi- 
modal distribution. The appearance of a larger value of 
r ^ 3 X 10^ s at 538 K reflects the presence of the notice¬ 
ably slower relaxation connected with the aging process. 

Since XPCS measures the contribution from all com¬ 
ponents of the system, according to their characteristic 
time scale of motion, it might be speculated if this behav¬ 
ior is attributed to the decoupling of individual compo¬ 
nent diffusivities, as observed in the supercooled liquid by 
Bartsch et al [21]. We rule out this hypothesis, as those 
self-diffusion data could lead only to a collective dynam¬ 
ics given by the average of distinct correlation functions 
with decay times very far from one another, and not to a 
temporal intermittent dynamics. The avalanche-like and 
intermittent temporal fluctuations in the collective dy¬ 
namics of Fig. emanate instead from distinct changes 
in the decay of correlations at the atomic length scale and 
thus contradict our intuition of a steady slowing down 
of the glass’s internal clock during aging. Intermittent 
dynamics, however, was not observed in earlier XPCS 
measurements of a Mg 65 Cu 25 Yio MG following a similar 
thermal protocol [9]. This suggests a possible connec¬ 
tion to the higher fragility of the Pd 43 Cu 27 NiioP 2 o super¬ 
cooled liquid [22] , although further experiments on other 
MG systems would be necessary to establish a definitive 
relationship. 

In Fig. [^the evolution of the TTCF is monitored af¬ 
ter heating the now partially aged sample directly to 563 
K [14]. A long regime of fast, stationary dynamics is 
observed up = 5 x 10^ s, which is interrupted by a 
sudden, then gradual broadening of the TTCF up to 
tw = 1-1 X 10^ s [Fig. Qa)]. The ^ 2 (^ 0 ,^) calculated 
in the initial stationary regime displays a fast dynamics 
with r = 150 s and a compressed exponential shape pa¬ 
rameter oi (3 = 1A [Fig. Qb)]. Afterwards, the dynamics 
has slowed considerably as the sample ages (r = 960 s), 
accompanied by a drop in the value of (3 to 0.9. Dur¬ 
ing longer annealing we observe the dynamics to speed 
up (r = 500 s) and the shape of the correlation function 
becomes re-compressed with p = 1.2 [Fig. Bc)^§d)]. 

The unique compressed shape of the correlation func¬ 
tions has been linked to ballistic-type motions that arise 
from the presence of internal stresses stored into glassy 
(or jammed) states during structural arrest [2l [9l [T^ [23] . 
This is in contrast to diffusive processes, characteristic of 
mass transport in liquids, which are driven by thermal 
energy and result in simple or stretched (/3 < 1) expo¬ 
nential decays [24] . By choosing (3 as an ad-hoc metric to 
reflect the evolution of microscopic stresses in the sample 
during aging, we can thus relate the sudden increase in r 
and concomitant decrease in P observed in these measure- 
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FIG. 4. (a) Measured TTCF after directly heating the annealed sample to T = 563 K. (b) Normalized one-time correlation 
functions calculated in the initially stationary (hlled squares) and aging (hlled circles) regime, (c) TTCF corresponding to the 
continuation of the measurement at T = 563 K. (d) Evolution of r and [3 over the entire length of the measurement. 


ments (see also El) to transitions from more-stressed to 
less-stressed states as the sample ages [2]. In fact, the 
decrease in P from 1.5 to ^ 1.2 [inset, Fig. ib)] agrees 
very well with the short- and the long-time values, respec¬ 
tively, predicted by the model of Bonchand and Pitard 
explaining the compressed decay of the correlation func¬ 
tions in terms of a release of internal stresses [23] . In this 
sense, we envision aging in MGs as the consequence of 
the gradual release of atomic-level stresses. 

A structural mechanism has been proposed that links 
these rearrangements to localized micro-collapses of 
groups of particles, which trigger subsequent collapses 
in neighboring regions through the formation of stress 
dipoles [23|. This avalanche-like dynamics has been 
found in a number of phenomena including martensitic 
transformations m, deformation of MGs [20l|26], crys¬ 
tallization of a hard-sphere glass m and shear flow of 
droplet emulsions through a thin opening [28] . We there¬ 
fore believe that these microscopic rearrangement events 
act as important mediators in the aging of MGs via the 
cooperative relief of atomic-level stresses. This point is 
illustrated in Fig. [^ where our heating of the annealed 
sample to 563 K apparently reduces the degree of in¬ 
ternal stresses, allowing the aging process to eventually 
re-activate. 

The transition in Fig. IT between stationary and ag¬ 


ing dynamics is accompanied by a decrease in (3 from 1.4 
to 0.9, thus suggesting a transition from ballistic to dif¬ 
fusive, liquid-like particle motion. Additionally, the final 
value of (3 = 1.2 measured at 563 K indicates that the 
sample is still stuck in a non-equilibrium state. This is 
surprising, as the duration of the annealing performed at 
this temperature 3 x 10^ s) far exceeds the extrapo¬ 
lated macroscopic r of the equilibrium liquid 1 x 10^ 
s). It would thus appear that - at least in the micro¬ 
scopic sense - full re-equilibration into the supercooled 
liquid is astonishingly difficult. This is further confirma¬ 
tion that macroscopic and microscopic aging in this MG 
clearly obey different laws. 

In closing, we note that similar aging dynamics has 
also been observed in colloidal gels and associated with 
intermittent ballistic particle motion [29]. Our results 
thus strengthen the similarities between soft matter and 
MGs, as well as their differences with respect to network 
glasses, which are characterized by a stretched exponen¬ 
tial decay of the density correlation function, even in the 
glass m- The intermittent dynamics reported here for 
the first time in a structural glass is indicative of a com¬ 
plex energy landscape m, characterized by numerous 
minima, in which the system dwells for long periods of 
time in stable configurations and reflect the presence of 
both localized and cascade relaxation dynamics, as in 
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the simulations of Ref. [20] . In this Letter we have also 
shown that, during the long annealing at each tempera¬ 
ture (which was not simulated in Ref. [20]), we observe 
a spontaneous re-activation of the aging from a local¬ 
ized dynamical regime, not reported in previous works. 
Nevertheless, these results should provide the impetus 
for a deeper exploration of the connection between local 
features of the PEL and out-of-equilibrium dynamics by 
means of simulation efforts, as well as for the develop¬ 
ment of new ways to exploit additional information from 
XPCS experiments. 

Karim L'Hoste and Hugo Vitoux are gratefully ac¬ 
knowledged for their technical support during the XPCS 
measurements. It is a pleasure to thank Yuriy Chuskin 
for providing the code for the XPCS data treatment. 
Z.E. is supported by the German Academic Exchange 
Service (DAAD). E.P. is supported by MINECO, grant 
EIS2014-54734-P and Generalitat de Catalunya, grant 
2014SGR00581. 


* Zachary.Evenson@frm2.tum.de 
[1] L. C. E. Struik, Physical aging in amorphous polymers 



plied Physics Letters 91, 051913 (2007): M. E. Launey, 
R. Busch, and J. J. Kruzic, Acta Materialia 56, 500 



825 (1959): A. Madsen, A. Fluerasu, and B. Ruta, 
in Synchrotron Light Sources and Free-Electron Lasers 
(Springer, 2014) pp. 1-21. 

[7] W. Gotze, Journal of Physics: Condensed Matter 11 , A1 
(1999) 

[8] W. Kob and H. C. Andersen, Physical Review E 52, 4134 
(1995) 

[9] B. Ruta, Y. Chushkin, G. Monaco, L. Cipelletti, 
E. Pineda, P. Bruna, V. M. Giordano, and M. Gonzalez- 
Silveira, Physical Review Letters 109, 165701 (2012) 

[10] M. Leitner, B. Sepiol, L. M. Stadler, B. Pfau, and 
G. Vogl, Nature materials 8, 717 (2009) 


[11] M. Leitner, B. Sepiol, L. M. Stadler, and B. Pfau, Phys¬ 
ical Review B 86, 064202 (2012) 

[12] B. Ruta, G. Baldi, G. Monaco, and Y. Chushkin, The 
Journal of Chemical Physics 138 , 054508 (2013) 

[13] B. Ruta, G. Baldi, Y. Chushkin, B. Ruffle, L. Cristo- 
folini, A. Fontana, M. Zanatta, and F. Nazzani, Nature 
communications 5, 3939 (2014) 

[14] “See the Supplementary Material for more information 
on the DMA measurements, data analysis and structural 
characterization.”. 

[15] C. A. Angell, Science 267 , 1924 (1995) 

[16] G. J. Fan, H.-J. Fecht, and E. J. Lavernia, Applied 
Physics Letters 84, 487 (2004) 

[17] C. A. Angell, K. L. Ngai, G. B. McKenna, P. F. McMil¬ 
lan, and S. W. Martin, Journal of Applied Physics 88, 
3113 (2000) 

[18] 1. M. Hodge, Journal of Non-Crystalline Solids 169, 211 
(1994) 

[19] B. Doliwa and A. Heuer, Physical Review Letters 91 , 
235501 (2003), arXiv:0306343 [cond-mat]: A. Heuer, 
Journal of physics. Condensed matter : an Institute of 
Physics journal 20 , 373101 (2008) 

[20] Y. Fan, T. Iwashita, and T. Egami, Physical Review 
Letters 115 , 045501 (2015) 

[21] A. Bartsch, K. Ratzke, A. Meyer, and F. Faupel, Phys¬ 
ical Review Letters 104 , 195901 (2010) 

[22] J. C. Mauro and R. J. Loucks, Physical Review E 78, 
021502 (2008): S. Wei, Z. Evenson, 1. Gallino, and 
R. Busch, Intermetallics 55 , 138 (2014) 

[23] J. P. Bouchaud and E. Pitard, European Physical Journal 
E 6, 231 (2001), arXiv:0107465 [cond-mat]: E. E. Ferrero, 
K. Martens, and J.-L. Barrat, Physical Review Letters 
113 , 248301 (2014), arXiv:1407.8014 

[24] A. Meyer, Physical Review B 81 , 012102 (2010): F. Yang, 
T. Kordel, D. Holland-Moritz, T. Unruh, and A. Meyer, 
Journal of Physics: Condensed Matter 23 , 254207 (2011) 

[25] L. Muller, M. Waldorf, C. Gutt, G. Griibel, A. Mad¬ 
sen, T. R. Finlayson, and U. Klemradt, Physical Review 
Letters 107 , 105701 (2011); C. Sanborn, K. F. Ludwig, 
M. C. Rogers, and M. Sutton, Physical Review Letters 
107 , 015702 (2011) 

[26] J. Antonaglia, W. J. Wright, X. Gu, R. R. Byer, T. C. 
Hufnagel, M. LeBlanc, J. T. Uhl, and K. A. Dah- 
men. Physical Review Letters 112 , 155501 (2014): J.- 
O. Krisponeit, S. Pitikaris, K. E. Avila, S. Kiichemann, 
A. Kruger, and K. Samwer, Nature communications 5 , 
3616 (2014) 

[27] E. Sanz, C. Valeriani, E. Zaccarelli, W. C. K. Poon, M. E. 
Cates, and P. N. Pusey, Proceedings of the National 
Academy of Sciences of the United States of America 
111 , 75 (2014) 

[28] D. Chen, K. W. Desmond, and E. R. Weeks, Soft Matter 
8, 10486 (2011), arXiv:1105.3099 

[29] A. Duri and L. Cipelletti, Europhysics Letters 76 (2006), 
10.1209/epl/i2006-10357-4, arXiv:0606051 [cond-mat] 







